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ABSTRACT 

We use several main-sequence models to derive distances (and extinctions), with sta- 
tistically meaningful uncertainties for 11 star- forming-regions and young clusters. The 
model dependency is shown to be small, allowing us to adopt the distances derived 
using one model. Using t hese d istances we have revised the age order for some of the 
clusters of iMavne et al.l (|2QQ7l ). The new nominal ages are: ~ 2 Myrs for NGC6530 
and the ONC, w 3 Myrs for A Orionis, NGC2264 and a Orionis, ^4-5 Myrs for 
NGC2362, « 13 Myrs for h and x Per, « 20 Myrs for NGC1960 and « 40 Myrs for 
NGC2547. In cases of significantly v ariable extinction we hav e derived individual ex- 
tinctions using a revised Q-method ([Johnson fc Morganlll953l ). These new data show 
that the largest remaining uncertainty in deriving an age ordering (and necessarily 
ages) is metallicity. We also discuss the use of a feature we term the R-C gap overlap 
to provide a diagnostic of isochronal age spreads or varying accretion histories within 
a given star-formation-region. Finally, recent derivations of the distance to the ONC 
lie in two groups. Our new more precise distance of 391^g 2 pc allows us to decisively 
reject the further distance, we adopt 400 pc as a convenient value. 

Key words: starsievolution - starsiformation - stars: pre-main-sequence - tech- 
niques: photometric - catalogues - (stars) Hertzsprung-Russell H-R diagram 



1 INTRODUCTION 

Colour- magnitude diagrams (CMDs) of star- format ion re- 
gions (SFRs) provide, in combination with model isochrones, 
an excellent tool with which to determine distances, ages 
and individual stellar masses. These parameters are criti- 
cal for determining initial mass functions (IMFs) for stellar 
populations and discovering the possible impacts of local en- 
vironment (such as the effect of ionising winds from massive 
stars) on disc lifetimes and on star and planet formation and 
evolution. Many calculations of IMFs, disc fractions etc are 
available but they are derived in heterogeneous ways, thus 
hints of the effects of environment are only recently begin- 
ning to emerge (e.g. iMavne et all 120071 : IStolte et al1l20Q4l ). 

Very precise photometry (« 1%) is routinely avail- 
able, along with sophisticated stellar models. However, cur- 
rent parameter derivations from CMDs still have relatively 
large uncertai nties and are model dependent (see the dis- 
cussions in iBonatto et all 120041: IPinsonneault et all 120041 : 
iNavlor et all 120021 : IMavne et all l200/t ). Thus current age 
(and distance) uncertainties all but 'wash-out' any envi- 
ronmental effects. Clearly more robust constraints would be 
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available for current stellar theories if more precise parame- 
ters could be extracted from the CMDs of SFRs. 

In IMavne et al . (2007) we created an age ladder for a 
range of pre-MS populations. The first stage was to create 
empirical isochrones by fitting splines to the pre-MS locus. 
Overlaying them in absolute magnitude and intrinsic colour 
results in an age ladder, with the youngest SFRs at the 
brightest absolute magnitudes. SFRs with almost indistin- 
guishable positions in the CMD were grouped, and nominal 
ages assigned to each group. Thus the age sequence (though 
not the nominal ages) are free from the proble ms associated 
with p re-main-sequence (pre-MS) models. In IMavne et al.l 
(2007) we had to adopt literature distances for the studied 
SFRs. These distances were derived using a range of differ- 
ent methods and their uncertainties proved to be the largest 
remaining contributor to the uncertainties in our age ladder 
placements. Previous distances have chiefly been derived us- 
ing main-sequence (MS) isochrone fitting, pre-MS isochrone 
fitting or from HIPPARCOS parallax measurements. MS 
isochrone fitting provides distances based on the positions 
of MS stars in a CMD, which are independent of uncertain- 
ties in age. Pre-MS isochrone fitting also uses the positions 
of stars in a CMD, but in this met hod the derived distance s 
are degenerate with age (see e.g. INavlor &; Jeffries! l2006h . 
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Finally, distances derived from HIPPARCOS parallax mea- 
surements are only available for a few SFRs included in this 
paper, a Ori, NGC2547 and A Ori, with all except A Ori 
having large uncertainties. 

Of those methods used to derive distances, the most 
suitable for the derivation of an age ladder is clearly MS 
isochrone fitting. Fully convective pre-MS stars (in young 
SFRs) are separated in a CMD from those stars on the MS 
which have radiative cores. The transistion region or gap 
in the CMD (measurable in c olour), we term the radiative- 
convective gap (R-C gap, see iMavne et al.l 120071 for intro- 
duction of term and discussion). Once stars have crossed the 
R-C gap their position in a CMD is almost independent of 
age until they reach the turn off. However MS fitting has not 
yielded the precision one would expect in distance estimates. 
This is due to two significant problems. Firstly the position 
of MS isochrones, although temporally static, is model de- 
pendent, with different studies adopting different models. 
Secondly distances are most often derived using 'by eye' fit- 
ting of models to the data, yielding ill-defined uncertainties. 
A full discussion of previo us fitting methods can be found 
in lNavlor fe Jeffries! jiooih . 

In this paper we solve both of these problems. Firstly 
we show that the model dependency is small for the model 
isochrones studied. We then adopt the distances from a sin- 
gle MS model. This allows us to derive a set of precise 
distances, accurate relative to each other, which we term 
"relative distances". Sec ond we use the r 2 fitting technique 
(|Navlor fc Jeffr ies 2006), a new rigorous and self-consistent 
method of fitting stars to isochrones, which yields statis- 
tically meaningful uncertainties. This presents us with the 
opportunity to achieve more precise distances from the fit- 
ting of high-mass (HM) or MS stars. 

The rest of this paper is laid out as follows. In Section [2] 
we detail the literature sources, the nature of the data used 
and any sequence selection carried out on the stars. Sec- 
tion [3] details the different model isochrones and photomet- 
ric calibrations used. Section 2] describes the fitting process. 
This is done primarily by way of an example in Section 14.11 
Section [5] describes the derivation of individual extinctions, 
in particular Section [5?2l describes an revised Q-method for 
calculating approximate individual extinctions. The results 
for all the isochrone calibrations and methods are presented 
in Section [6] Section [7] outlines our results for one adopted 
model with our best-fitting distances given in Table [8] In 
section [8] we discuss the implications of the individual dis- 
tances to several key SFRs (Section 18. 1 p and those of the 
entire dataset. The implications of the dataset on metal- 
licity (Section I8.2.1[) , age spreads and the R-C gap overlap 
(Section I8.2.2[) and secular evolution within the SFRs with 
particular reference to disc fractions (Section I8.2.3|) are dis- 
cussed. The reader interested in distance and reddening val- 
ues should skip to Section [6] for the values derived from all 
the models used. 



2 THE DATA 

All the datasets presented in this work are from literature 
sources. To avoid distance-age degeneracy problems and to 
minimise the effect of age assumption on our distance deriva- 
tions we have only fitted HM stars on or near the MS. We 



have used the memberships adopted in the original source 
and made further photometric cuts. The sources of photo- 
metric data and initial memberships are shown in Table [1] 
The photometry is in the Johnson-Cousins system unless 
otherwise stated. 

Further photometric cuts are required to remove non- 
member stars which may have satisfied the membership cri- 
teria within the original publication. In addition photomet- 
ric cuts are required to select the correct part of the se- 
quence. Motion in a CMD as a function of age is rapid in 
the pre-MS and post-MS phases, compared to that on the 
MS. Therefore inclusion of stars in these phases would intro- 
duce an age dependency into our distances. Thus we make 
the two photometric cuts detailed below, which in all cases 
results in a clearly identifiable MS which lies clear of the 
contamination. 



2.1 The turn-off cut 

If we include stars which are too bright they may have 
evolved away from the MS (turn-off), so we make a photo- 
metric cut at the turn-off, for the nominal age, despite the 
fact that in some cases the MS appears to extend above the 
turn-off (see Section [8] for a discussion of this discrepancy) . 

There is also a shift in the position of the MS as the 
isochronal age increases, as the stars are 'preparing' to turn- 
off. We have examined this effect and the result of assuming 
the wrong age for the model sequence. The effect is negligible 
in our experiment as explained in Section T4. 2 1 

2.2 The turn-on cut 

If we include stars which are too red they may be pre-MS 
stars where age is degenerate with distance. The positions 
in colour at which stars join the MS (turn-on) and hence 
where the turn-on cut should be made is predicted by pre- 
MS isochrones. F igure [T] shows the pre-MS isochrones of 
ISiess et al.l ([2000h for ages of 1, 3, 5, 15 and 40 Myrs (the 
latter being our oldest SFR). The positions in colour of the 
turn-ons are shown in Figure [T] However for the younger 
SRFs the observed MS in a CMD often appears to extend 
redder and fainter than predicted by the pre-MS isochrones. 
So we have decided on the positions of the cuts empirically, 
that is to say we have identified the bottom of the MS in 
the data and placed the cut there, at the blue edge of the 
radiative-convective gap. Table [1] sho ws the posi t ions o f the 
cuts predicted from the isochrones of ISiess et al.l ((2000), the 
reddening (see S ect ions l4.1.2l and 171), the actual cut employed 
and the age assigned in IMavne et al.l ([20071 ) . Using this em- 
prical cut, as opposed to that from theory results in signifi- 
cantly more precise distances for some of our youngest SFRs. 
This is because the distance is primarily derived from the 
curve of the MS towards the red at fainter magnitudes. A 
discussion of the implications of the MS extending below its 
theoretical terminus can be found in Section [8] 



3 THE MODELS 

In this work we hav e used the MS s tellar interior 
models of the Padova (jGirardi et al.l l2002h and Geneva 



Fitting the young main- sequence 3 



SFR 


B — V 
Cut 


E(B - V) 


(B - 
Data 


V) Cut 
Theory 


Nominal age (Myr) 
Mavne et al. (2007") 


Data 
type 


Source 


Size 
(arcmins) 


the ONCW 






0.1 


-0.13 


1 


VI & T eff 


Hillenbrand (1997) 


20 


NPPR^n( 2 ) 


u. o 


n ^9 
u.oz 


u. io 


n 17 

-U. 1 l 




TT RV T 
U D V I 


ouiir et ai. { zuuu j 


io 


NGC2244( 2 ) 


0.45 


0.47 


-0.02 


-0.13 


^2( 3 ) 


UBVIH a 


Park & Sune (2002^ 


15 


NGC2264 


-0.02 


0.03 


-0.05 


-0.1 


3 


UBVRI 


Mendoza V. & Gomez (1980) 


30 


NGC2362 


0.04 


0.09 


-0.05 


-0.1 


3 


UBV 


Johnson & Morgan ( 1953) 


30 


A Ori 


0.20 


0.11 


0.09 


-0.1 


3 


UBV 


Murdin fe Penston (1977) 


30 


a Ori 


0.03 


0.06 


-0.03 


-0.06 


4-5 


V T B T ^ 


Caballero (2007) 


30 


X Per( 7 ) 


0.7 


0.50 


0.2 


0.2 


13 


UBV 


Slesnick et al. (2002) 


2 


h Per( ? ) 


0.74 


0.54 


0.2 


0.2 


13 


UBV 


Slesnick et al. C2002) 


2 


NGC1960 


0.25 


0.20 


0.05 


0.2 


16(5) 


UBV 


Johnson & Morgan ( 1953) 


10 


NGC2547 


0.1 


0.04 


0.06 


> 0.2 


38( 6 ) 


UBV 


Claria 1982) 


15 



Table 1. The empirical and theoretical B — V cuts for each SFR. The apparent cut has been conve rted into intrinsic colour using the 
E(B — V) derived in Section [7] Also shown is the nominal age for each SFR from Ma vne et al.l ( 2007), the data type used and source for 
the data and initial memberships. The notes are as follows. (1) Individual extinctions derived from T e ff, see Section [5] (2) Individual 
extinctions derived us ing the Q-me thod (see Section [5?2)) (3) Age from IPark & Sund (|2Q02h . (4) Data in the TYCHO photometric 
system. (5) Age from ISanner et aTl (|2000T ). (6) Age from iNavlor fc J effries (20( 3). (X) F or h and x P er stars have additionally been 
selected using apparent distance from the cluster centre fas defined in lMavne et aTl l2007). selecting a circular area around the cluster 
centre of 2 arcmins in radius. 
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Figure 1. 1, 3, 5, 15 and 20 Myr isochrones from ISiess et all 
( 2000) showing the MS, and pre-MS. B — V colour cuts are shown 
for each age to isolate the MS. These are: -0.17, -0.1, -0.06 and 
0.2 for 1, 3, 5, and older sequences respectively. 



([Leieune &; Schaererll200lh groups. These provide an effec- 
tive temperature (T e //), luminosity and surface gravity. 
These values must then be converted into colours and magni- 
tudes in the required photometric system (Johnson-Cousins) 
to allow the fitting of photometric data. Colours are found 
using a T e // to colour relation, and magnitudes using the 
bolometric correction to the luminosity. Both the colour- 
T e f f relation and bolometric correction come from using the 
parameters from a stellar interior model to find the correct 
model atmosphere and then folding the resulting flux dis- 
tribution through appropriate photometric filter responses. 
Once this is achieved the photometric colours and magni- 
tudes must then be calibrated to a standard scale, using the 
colours of Vega in the photometric system. We have used 
three main isochrone and extinction systems calibrated to 
two different Vega colour systems, and we now detail each. 



3.1 Geneva 

The Geneva isochrones (as provided in iLeieune &; Schaererl 
2001) are from the Geneva stellar interior models (basic set) 
in conjunctio n with the updated BaSeL-2.2 model atmo- 
spheres from IWestera et aD (|l999h . To derive photometric 
magnitudes (V) they hav e adop ted the bolometric correc- 
tions from ILeieune et al.l (Il998h w hich are defined to fit 
the empirical scale of iFlowerl (|1996l ) (not calibrated to the 
Sun). They calculate colours for Johnson- Cousins photom- 
etry using the filter response func tions of iBuser &; Kuruczl 
(| 19781 ) (UBV) and iBesselll (|l979h (RI). The colours and 
magnitudes of the isochrone are then calibrated to Vega 
colours of zero. For these isochrones we use the canonical 
extinction vectors, namely E(U - B)/E(B — V) = 0.73, 
A{V)/E{B -V) = 3.1 and A{V)/E{V - I) = 0.41. 

3.2 Geneva-Bessell 

For the Geneva-Bessell isochrones we h ave used the inte- 
rior models of ILeieune specifically their 
basic model set ("c") generally applicable for stars with 
M < 1 2M&. Their convers i on to photometric co lours follows 
that of lBessell et~ai1 (|l998h . lBessell et al.1 (Il998h u se the AT- 
LAS9 atmosphere models of ICastelli et al.l (Il997h fat solar 
metallicity only) and the filter responses of IBesselll (|l990h 
(UBVRI), to calculate the colour- T e // relations and bolo- 
metric corrections. The resulting Johnson-Cousins photom- 
etry is then calibrated to Vega colours of zero. In practice 
we carry out these conversions in our own code, since that 
allows us to also calibrate to what we term the non-zero 
system where (B — V)v ega — —0.002 and (U — B)vega 
().()() 1. 



3.3 Padova-Bessell 

For the P adova-Bessell i sochro nes we use the stellar interior 
models of lGirardi et al.l (|2QQ2h . The colours an d magnitudes 
are th en calculated using the conversions of iBessell et al.l 
(1998) as for the Geneva-Bessell isochrones. The resulting 
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Johnson- Cousins photometry is then calibrated to either 
Vega colours of zero or the non-zero system. 

3.4 Extinction vectors in the Bessell system 

It is well known that extinction vectors ar e actually a 
funct ion of intrinsic colour (or spectral type) . iBessell et al.l 
(1998s) provide extinction vector s as a functio n of colour 
based on the extinction curves of lMathisI (|l990h . We there- 
fore use these extinction vectors for the Geneva-Bessell and 
Padova- Bessell models. 

The extinctions are provided for an E(B — V) — 0.3. 
To check the range of extinctions to the SFRs studied here 
does not have a significant effect, we folded the solar abun- 
dance AT LAS 9 spectra with the " new" opacity distribution 
fu nction (iCastelli &; Kuruczl 120041 ) through the bandpasses 
of lBessell et al.l (|l998l ). We t hen reddened t he sp ectra ac- 
cording to the prescription of ICardelli et al.l (|l989h to yield 
an E(B — V) of ap proximately 0.3. (This funct ion is the one 
tabulated in the Mathis (1990) paper used by Be ssell et al.l 
(1998Q). We used an R v of 3.2 since, when folded through 
the bandpasses of lBessel3 {1990) we fo und this gave the bes t 
match to the BV extinction vector of Bessel l et al.l ([1998). 
We calculated the difference /E(B - V) for values of 

E(B - V) of 1 (typical of the SFRs we have fitted) and 3 
(approximating to the highest reddening of any SFR fitted). 
For B — V < 1.5 we find the largest difference is -0.02 mags, 
which has a negligible impact on our fits. 

3.5 TYCHO photometry 

For a Orionis the data were taken in the TYCHO photomet- 
ric syst e m. W e have fitted these data using the conversion of 
Bessell (2000) transform our Geneva-Bessell isochrones into 
the TYCHO system. We have defined extinction vectors in 
the TYCHO photometric system, in a process similar to 
that described in Section 13.41 We set our zero-points by re- 
quiring we re produced the T e ff vs B — V relationship of 
Besselll dl990h and the B — V vs A(B - V) relationship of 
Besselll (J2000). This gave, for (B - V) T < 0.065 



(A v )t 



-- 3.358 + 0.237(B - V>, 



E{B-V) T 
and for 0.065 < (B — V) T < 0.5 



3.387 - 0.207(B -V) T . 



(1) 



(2) 



E(B-V)t 

These fits never deviate from the calculated curve by more 
than about 0.01 mags. At B — V=0 the BV curves match 
to within 0.01 mags, though this worsens to 0.05 mags at 
B - V=l. 



4 THE FITTING METHOD 

Throughout this example and later sections all the isochrone 
fits displayed and tests undertaken use the Geneva-Bessell 
Vega-zero isochrones. These isochrones are also adopted in 
Section where we draw implications from the resulting 
distances and age ordering. As shown in Section [6] the model 
dependency between the different isochrones is practically 
very small. 



4.1 Example fit: x P er 

It is most instructive to describe the fitting method via 
an example. Here we us e the cluster v Pe r which is ap- 
proximately 13 Mv rs old (iMavne et alJl2007h . at a distance 
modulus of 11.85 dsiesnick et al.l l2QQ2h with an A v w 1.6 
(Ma vne et al.l 120071 ) (the extinction is reasonably uniform). 
First, we describe the derivation of distance. Deriving a dis- 
tance does require a known extinction, the derivation of 
which is described later in this section. 



4-1.1 Fitting statistic, distance derivation 

The fitti ng statistic used in this work is r 2 , which is intro- 
duced in iNavlor fe Jeffries! <l2006h. Fitting of MS d ata using 
this technique is described in Jeffries et al. (2007|). r 2 is es- 
sentially a generalised x 2 statistic including uncertainties in 
two dimensions, and models with a two-dimensional distri- 
bution as opposed to a single isochronal line. The best fitting 
model is found by minimising r 2 . 

Once we have selected an isochrone we use a Monte- 
Carlo method to generate a colour- magnitude probability 
grid. Each pixel in this grid is assigned a value that gives 
the probability of finding a star drawn from the popula- 
tion represented by the isochrone, at any given colour and 
magnitude. The model is then adjusted through a range 
of distances and the values of r 2 for each star summed to 
calculate a total value of r 2 for each distance step, as de- 
tailed in Naylor & Jeffries! (2006). These r 2 contributions 
were clipped, i.e. the contribution to the total r 2 for any 
single data point value is capped at some set value. Clipping 
avoids erroneously included non-members or anomalous ob- 
jects many a from a given model overwhelming the result. 
The lowest total r 2 was then selected as the best fitting 
model. 

Once a given fit was completed a probability of obtain- 
ing the resulting r 2 , P r (r 2 ), was calculated. If P r (r 2 ) is far 
from 50%, the fitting is repeated with an additional sys- 
tematic uncertainty added to the data. This is analogous to 
enlarging error bars to achieve a xl — 1- This process en- 
sures the model is a good fit to the data and we are then 
able to derive the parameter uncertainties. To derive these 
uncertainties, we used the bootstrap method described in 
INavlor &; Jeffries! (|2006l ), repeating the fitting 100 times and 
deriving 68% confidence intervals. 

The resulting fit for x P er is shown in Figure [2] The 
derived distance modulus and 68% confidence interval is 
11.79 < 11.83 < 11.88. This agrees wi th the most recent lit - 
erature derivation of 11. 85 ±0.05 from Slesnick et al.1 <|2002h , 
which was the same dataset. 

Figure [2] and all subsequent figures showing fitted data 
have several elements requiring explanation. The shaded 
area shows the probability density of finding a star at a 
particular colour and magnitude (or colour and colour for 
extinction fitting) . This density, p is that from Equation 1 of 
INavlor &; Jeffries! (|2006h generated for a specific isochrone, 
in this case Geneva-Bessell. The circles show the positions of 
the photometry and the bars give the uncertainties in mag- 
nitude and colour. For all the figures showing fitted data 
(except those using individual extinctions, see Section [5} 
the models have been adjusted to the natural space of the 
data i.e. apparent colour and magnitude. 
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Figure 2. The distance fit for x P er to the Geneva-Bessell 
isochrones. See Section PTTl for details of the symbols. 
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Figure 3. Figure showing E(B — V) fit of % Per to the Geneva- 
Bessell isochrones. See Section 14.1.21 for details of the symbols. 
Here a systematic shift in E(B — V) as a function of colour 
(U — B) is evident, as in Figure [9] T his is probably due to a 
difference in the photometric systems of Bessell et al.l (119981 ) and 
ISlesnick et"aH J2002) (see Section E2J). 



4-1.2 Mean extinction fitting 

To allow us to derive a distance, an extinction is required 
and is indeed crucial as changing it will change the distance 
derived for the stars we are fitting. We can derive a mean 
extinction by fitting the data to an isochrone in a colour- 
colour diagram. Where we have UBV photometry we have 
simply fitted the sequence in U — B vs B — V in a similar 
fashion as that for a distance. However, instead of chang- 
ing the distance, we evaluate r 2 at different values of the 
reddening. The resulting fit can be seen in Figure [3] with 
a best fitting E(B - V) = 0.50. Figure [3] and all the sub- 
sequent figures showing extinction fitting contain the same 
components as those for distance fitting. 

In general fitting for an extinction using the r 2 method 
is desirable as photometric individual extinction methods 
rely on the star being a sinlge star or an equal mass-mass 
binary (see in Section \5.2^ . However, in some cases the dis- 
persion in an E(B — V) fit is too large, i.e. the scatter 



around the isochrone in the U — B vs B — V is too large 
to confidently assign one mean extinction. This is the case 
where there is significantly variable reddening across a SFR. 
Here we are forced to derive reddenings and therefore extinc- 
tions for each star from UBV photo metry alone using the 
Q-method ([ Johnson &; Morgan! 1 19531 ). This is the case for 
NGC6530, NGC2244 and A Orionis. In addition in the case 
of the ONC UBV colours are not available, here alterna- 
tive methods must be used to derive extinctions. All these 
individual extinction derivations are detailed in Section [5] 

4.2 Practical effect of assumptions 

During fitting we have made three main assumptions about 
the MS and it is important to examine the validity of these. 
Firstly we have used an isochrone at an assumed nominal 
age which could be inaccurate for a given SFR. This may 
be important both just prior to turn-off and at the turn- 
on. Secondly, by using MS isochrones we have implicitly as- 
sumed that none of the stars fitted are on the pre-MS. Lastly, 
we have used isochrones of approximately solar composition 
(Z = 0.02). 

4.2.1 Age assumption 

As a (coeval) stellar population ages, stars of a decreas- 
ing mass turn-off from the MS. In addition, stars close to 
but nominally below the turn-off age also move slightly red- 
wards in position in a CMD. As we remove stars brighter 
than the turn-off fSection !2.ip we avoid the age dependency 
of this feature, but we still retain the age dependency of the 
upper-MS. This shift of the upper-MS is small and further- 
more in principle is modeled as we use a MS at the nominal 
age. However it is important to quantify this effect to ensure 
our distances are robust against assuming an incorrect age. 

To test this upper-MS age dependency we have simu- 
lated photometry of a population based on a lMyr Geneva- 
Bessell isochrone. We have then fitted these data (as de- 
scribed in Section UT} across a typical colour range, specifi- 
cally blue- ward of (B — V)o = 0.2 (as shown in Table [T]cuts 
in all our SFRs are blue- ward of this), to a lOMyr Geneva- 
Bessell isochrone. This provides a measure of any extra un- 
certainty one accrues if the age assumption is incorrect. The 
resulting distance modulus and reddening are dm ~ 0.001 
and E(B - V) « -0.028, for this factor ten in age. The ef- 
fect of assuming an incorrect age on a derived distance and 
E(B — V) is therefore negligible even for a large error in 
assumed age. 

4.2.2 MS isochrones for a possible PMS population 

As pre-MS stars approach the MS, just prior to the onset of 
hydrogen burning they enter a quasi-equilibrium state as hy- 
drogen ignition starts. This delays their arrival onto the MS. 
The rate at which this phase progresses is a function of stel- 
lar mass and is, over the mass range of interest, very short 
in comparison to the evolution of the pre-MS. However, as 
such stars are slightly brighter than MS magnitudes, the dis- 
tances derived may be systematically reduced as a function 
of age. 

Some stars in our selected sequence may still be in this 
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pre-MS phase even though appearing to be on the MS. This 
effect is greatest for the younger SFRs. However within the 
colour cuts we have used the deviation from the MS is a 
maximum of ~ 0.05 mag (in magnitude) and only affects a 
limited number of stars. Thus its overall effect will be much 
smaller than 0.05 mags. 

4-2.3 Composition 

As compositions are not available for all the SFRs stud- 
ied we have assumed solar metallicity (Z = 0.02). However 
there is evidence to suggest that some SFRs have a metallic- 
ity as low as half-solar (see Sect ion f8. 2.1)) . We must therefore 
quantify the possible effect varying composition could have 
on our relative distances. To do this we have again simulated 
a population, at lOMyr, using the Geneva isochrones with 
Z = 0.008 (the closest match in the library to half-solar) 
and fitted these data for extinction and distance with solar 
composition isochrones (as described in Section H~T]) . The re- 
sulting differences are, for reddening E(B — V)— 0.005 and 
for distance modulus, dm w 0.41. This is a negligible differ- 
ence in reddening and therefore extinction but a significant 
error in distance modulus. We have also fitted % Per using 
the half-solar metallicity isochrones. Here we must use the 
colours and photometric system provided with the Geneva 
isochrones as our Geneva- Bessell isochrones utilise atmo- 
spheres of solar metallicity. The resulting distance modulus 
is 11.32 < 11.36 < 11.41 and a reddening of E(B-V) = 0.51 
(using uncertainties of 0.018 for a P r (r 2 ) « 0.5). This means 
that if the metallicity is indeed Z ~ 0.01 the distance mod- 
ulus to x P er is actually ^0.5 mags less, making the stars 
older. This could have a major effect on the distances and 
ages of star-forming regions. Therefore composition informa- 
tion is vital in the future for more accurate SFR parameters. 
A further discussion of this problem can be found in Section 

m 



5 INDIVIDUAL EXTINCTIONS 

As stated in Section 14.1.21 it is sometimes not possible or 
desirable to derive a mean extinction using the r 2 fitting 
method. In this Section we detail the cases where individual 
extinctions have been derived. Source by source extinctions 
have been used for six SFRs, using two methods. In all of 
these cases the extinctions applied will change the relative 
positions of the stars in colour magnitude space, as opposed 
to applying a mean extinction where the entire dataset is 
simply translated (shifted). Therefore, in the case of individ- 
ual extinctions these must be applied prior to fitting, and 
the resulting figures show the data and model in intrinsic 
colour and extinction free magnitude. 

5.1 Effective temperatures 

5.1.1 the ONC 

For the ONC UBV photometry is unavailable, so fitting for a 
mean extinction using t he method detaile d in Section l4.1.2l is 
not possible. However iHillenbrandl (|l997T ) provides VI pho- 
tometry an of the centre of the ONC. We have used these 
temperatures, assumed a surface gravity (log(g)) of 4.5 and 



interpola ted to obtain intri nsic colours from the T e // rela- 
tions in lBessell et all (| 19981 ). This allows us to derive an Ay 
from E(V - I) using Ay = ((3.26 + 0.22(£ - V) )/(1.32 + 
0.06(y — I)p)) x E(V — I) (derived from the extinction vec- 
tors of Bessel fet al.lll998T ) and therefore calculate an unre- 
denned magnitude. These resulting magnitudes and colours 
can then be fitted to derive a distance as in Section 14.11 For 
the analysis in Section [7]of this paper we have fitted the data 
using the intrinsic B — V values appropriate for the effective 
temperatures, but we have also fitted the V — I data, the 
result of which can be found in Section [6] 

5.1.2 h and x Per 

The check the veracity of this technique, we compared the 
distances it yields for h and x P er 5 with t hose from Sec- 
tions |4] and [6] We used the T e / / data from Sles nick et al.l 
(2002) to derive intrinsic colours. However, as the T e // in- 
formation is only available for the hottest stars, the colour 
range for fitting is small. Practically this means the uncer- 
tainties are large, particularly towards greater distances, but 
the answers for both techniques are consistent. As a fur- 
ther test, we have supplemented the objects with T e // data, 
with stars from the main catalogue dereddened to intrinsic 
colours using the mean extinction from the r 2 fitting. All 
the two techniques give consistent distances (Table [5}, but 
those using the combined photometric and T e f / datasets are 
more precise (ranges of 0.01 and 0.04 mags for h and x P er 
respectively). This confirms that the T e // technique gives 
consistent answers, but we are wary of adopting the more 
precise technique for the remainder of this paper since our 
primary aim is to derive distances for many SFRs using a 
single method. 

5.2 Q-method 

The Q-method is a widely used method to derive individual 
source-by-source reddenings and therefor e extinctions using 
a U — B vs B — V colour- colour diagram. Ijohnson &; Morgan! 
(1953) model the MS as the straight line in Figure 2] As this 
Figure shows, this can result in errors in the derived Ay of 
up to 0.1 mag, to which must be added further 0.08 mag due 
to using colour independent extinction vectors. These are 
not intrinsic failings of the method, fitting a straight line to 
a modern isochrone, and using colour-dependent extinction 
vectors yields 

Q = 0.24(B - V) 2 + 3.257(5 - V) + 0.015, (3) 

in the Bessel fet al.l (Jl998) system. Within the ranges of 
colour in Table [2] this only differs by < 0.005 mag from 
interpolating onto the isochrone. However in our "revised 
Q-method" we derive reddenings using isochone interpola- 
tion. 

Implicit in using a model, are assumptions as to the age 
and metallicity of the isochrone. We have tested the effect of 
both of these, and provided one remains within the colour 
ranges given in Table [2] their effect on the derived Ay is 
less than 0.005 mags. The blue limit in Table [2] corresponds 
to the turn-off, the red limit to approximately U — B — 
B — V — 0. A further possible concern is shown in Figure 
[4] once the intrinsic colour is red ward of U — B ?^-0.2, there 
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Age 


U - 


B 




B 


- V 




1 


-0.17 < U - 


B < 


-1.15 


-0.06 < B 


-V < 


-0.30 


3 


-0.17 <U - 


B < 


-1.12 


-0.06 < B 


-V < 


-0.29 


5 


-0.17 <U - 


B < 


-1.08 


-0.06 < B 


-V < 


-0.28 


15 


-0.17 < U - 


B < 


-0.90 


-0.06 < B 


-V < 


-0.23 


30 


-0.17 < U - 


B < 


-0.72 


-0.06 < B 


-V < 


-0.22 



Table 2. The approximate validity range for Equation [3] (not 
including multiple valued solutions) for several ages derived from 
the Geneva-Bessell isochrones, denned as regions where the real 
isochrone moves more than 0.01 in colour away from the straight 
line model. 



is the possibility that the star actually lies on a redder part 
of the isochrone, leading to an ambiguity in the extinction. 
However, if the object is on the wrong part of the isochrone, 
the extinction is clearly anomalous, provided that the scatter 
in extinctions between stars is small. 

The major concern when calculating individual redden- 
ings and therefore extinctions using a colour-colour diagram 
is the effect of binarity. Figure U shows the range of colours 
occupied by unequal-mass binaries. In the absence of multi- 
plicity information, any colour-colour method must dereden 
stars onto the single star /equal- mass binary sequence (see 
Section I4T21) . As can be seen in Figure [3] without multi- 
plicity information it is unclear whether a star should lie on 
the single star /equal- mass binary sequence or in the region 
occupied by unequal-mass binaries. The effect and range of 
this problem is shown in Figure 2) The outer binary envelope 
has been modeled and the possible differences in derived Ay 
found. Figure [4] shows that binarity has an effect of up to 
AAy ~ 0.15 on the derived extinction. This effect becomes 
increasingly significant as one moves down (redder in U — B) 
the MS isochrone. 

As the revised Q-method does not account for the scat- 
ter binaries produce in a CMD, or colour-colour diagram, it 
is a statistically ill-defined process. Effectively most of the 
intrinsic scatter from the binary sequence and photomet- 
ric uncertainties is removed, in addition to that caused by 
variable extinction. Therefore, only when there is good ev- 
idence, from the mean extinction fitting method, that the 
extinction in a SFR is large and variable do we apply the 
revised Q-method. Thus we formulate the null hypothesis 
that the reddening or extinction is uniform. We fit to de- 
rive a mean extinction and subsequently a distance. These 
results can be found in Table [6] and are discussed in Section 
[6] However, in some cases of distance fitting, after applying 
a mean extinction, the addition of large systematic uncer- 
tainties was required to return a Pr(r 2 ) « 50%. In these 
cases we are forced to reject the null hypothesis and use the 
revised Q-method to derive individual extinctions. We be- 
lieve that additional systematic uncertainties of up to 2% are 
credible, therefore we apply the revised Q-method in cases 
where our added systematic uncertainties exceed this level. 

5.3 SFRs with significantly variable reddening 

In Figure [6] we show the r 2 mean extinction fit for NGC6530. 
There is clearly a large scatter, which is also reflected in the 
corresponding distance fit, Figure [5] This scatter improves 
significantly when we use the revised Q-method, as shown 
in Figure 1201 In addition the uncertainties in distance are 



0.10 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 

(B-V) 



Figure 4. Geneva-Bessell 10 Myr isochrone including binaries 
(gray region enclosed by black line) and the Q-method MS 
straight line (red line). The horizontal dashed lines are the points 
at which the Geneva-Bessell isochrones evolve away from the MS. 
The horizontal bold lines are the validity region of the original Q- 
method. The two angled dotted lines are the extin ction vectors 
at ea ch end of the Q-method validity range from iBessell et al.l 
( 1998). Finally the large dashed line shows the region below which 
the solution for a given star can be multi-valued. The inset shows 
the differences in derived extinction (A Ay)- The top dashed line 
shows AAy in the sense of the revised Q-method minus the old Q- 
method. The lower dashed shows AAy derived from the reddest 
binaries minus the old Q-method. The bold line is AAy derived 
from the reddest binaries minus the single star sequence. 

significantly smaller when using the revised Q-method. The 
same arguments apply for NGC2244 and A Ori, and in Sec- 
tion [7] all three clusters are fitted using extinctions from the 
revised Q-method. For completeness the parameters derived 
for these three clusters using both methods are presented 
in Section [6] The resulting distance moduli derived using 
a mean extinction or the revised Q-method are consistent 
within the uncertainties. 

We also attempted to use the revised Q-method for h 
Per, as it satisfied our criteria for non-uniform redenning. 
However, we found the E(B — V) to vary systematically 
as a function of colour. This systematic shift in extinction 
with colour is evident in Figure [9] It shows that to fit the 
hotter stars and cooler stars simultaneously would require a 
change in the gradient of the isochrone. The same trend is 
observed for \ P er m Figure [3] Therefore, we attribute this 
beha viour to difference s be tween the ph o tomet ric systems 
of ISlesnick et al.l (|2002h and lBessell et al.l (|l998h . Moreover, 
as this systematic shift in E(B — V) is not present in V vs 
B — V fit (see Figure [2} we further constrain the problem to 
a difference dominated by the U band. 



6 MODEL DEPENDENCY 

Our results for all models are given in Tables [3d The fits 
have been optimised by adjusting the systematic uncertain- 
ties in colour and magnitude such that the Pr(r) ~ 50% 
(actually 44%-66%). The values for the reddening, the dis- 
tance modulus with 68% confidence limits and, the added 
systematic uncertainties are shown. Table [3] shows the dis- 
tances derived for a Orionis after conversion to the TYCHO 
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Figure 5. NGC6530 distance fit using mean extinction, the fit 
for which is shown as Figure [6] 
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Figure 6. NGC6530 E(B - V) fit. 



photometric system of Bessell (2000). Tables |4] and [5] con- 
tain the results fitting using a direct conversion from T e f / to 
colour. The distances derived for each SFR using the mean 
extinction and revised Q-method for each set of isochrones 
are shown in Tables [6] and [7] 

These results allow us to perform a brief comparison 
of the isochrones we have used, although not the main aim 
of this paper, it may aid the reader in adopting a particu- 
lar result for a SFR of interest. Figure shows the derived 
distances and associated uncertainties (68% confidence in- 
tervals) for each SFR and each set of isochrones, where a 
mean extinction has been derived. Figure [7] shows that for 
any given SFR the scatter between models is smaller than 
the uncertainties from the data. However, one could argue 
for a systematic shift of « 0.05 in distance modulus depend- 
ing on the choice of Vega zero point. 



Geneva Bessell (0) 
Padova Bessell (0) 




Figure 7. The distances and uncertainties (68% confidence in- 
tervals) derived for all SFRs and all isochrones, where a mean 
extinction has been derived. The results shown for each cluster 
are for, from left to right, Geneva (0) (blue), Geneva Bessell (0) 
(black), Padova Bessell (0) (red), Geneva Bessell (/=0) (green) 
and Padova Bessell (/=0) (yellow). All the resulting distances 
agree within the uncertainties for all models. 



TYCHO photometry 


SFR 




a Ori 






Model (Vega cal) 




dm 


Adm 


Unc 


Geneva-Bessell(O) 




7.84 < 7.94 < 8.10 


0.26 


0.00 


Padova-Bessell(O) 




7.72 < 7.93 < 8.03 


0.31 


0.00 


Geneva-Bessell(/ = 


0) 


7.79 < 7.94 < 8.07 


0.28 


0.005 


Padova-Bessell(/ = 


0) 


7.72 < 7.93 < 7.98 


0.26 


0.005 



Table 3. The distance moduli derived for a Ori in the TYCHO 
photometric system for the different models, and calibrations to 
Vega colours of zero or Vega B-V = -0.002 and U — B = -0.004. 
Uncertainties required to achieve 0.40 < P r (r 2 ) < 0.60 are also 
shown. 



Intrinsic 


colours from T e f f 






SFR 


The ONC 




Model (colour used) 


dm 


Adm 


Unc 


Geneva-Bessell(B - V) 


7.91 < 7.96 < 8.03 


0.12 


0.010 


Geneva-Bessell(y — I) 


7.88 < 7.96 < 8.00 


0.12 


0.013 



Table 4. The distance moduli derived for the ONC using in- 
trinsic colours derived from T e ff values of iHillenbrandl (|l997T ) in 
different colour indices. Uncertainties required to achieve 0.40 < 
P r (r 2 ) < 0.60 are also shown. 
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Intrinsic colours from T e 


//' B ~ 


V 






SFR 

Model (stars used) 


X Per 

dm 


Adm 


Unc 


h Per 

dm 


Adm 


Unc 


Geneva-Bessell(Hot) 
Geneva-Bessell(Hot and cool) 


11.70 < 11.78 < 11.93 
11.83 < 11.84 < 11.87 


0.23 
0.04 


0.003 
0.00 


11.76 < 11.78 < 11.94 
11.83 < 11.83 < 11.84 


0.18 
0.01 


0.004 
0.019 



Table 5. The distance moduli derived for h and x P er using intrinsic colours derived from T e ff values of ISlesnick et al . (2002) alone, 
or supplemented by cooler stars dereddened using an average extinction. Uncertainties required to achieve 0.40 < P r (r 2 ) < 0.60 are also 
shown. 



Reddening from U BV fitting 



O 



SFR 




NGC6530 






NGC2244 






NGC2264 








Model (Vega cal) 




dm 


Adm 


Unc 


&(B-V) 


dm 


Adm 


T T„ _ 

Unc 


&(B-V) 


dm 


Adm 


T T„ _ 

Unc 


TP 

U(B-V) 




Geneva(O) 




10.29 < 10.38 < 10.48 


0.19 


0.0330 


0.35 


10.62 < 10.78 < 10.94 


0.32 


0.0470 


0.47 


9.15 < 9.40 < 9.53 


0.38 


0.018 


0.06 




vjrciic v a JJcoisciii u J 




10.16 < 10.34 < 10.43 


0.27 




0.32 


10.68 < 10.89 < 11.09 


0.41 


0.0450 


0.46 


9.26 < 9.37 < 9.52 


0.26 


0.0160 


0.04 




Padova-Bessell(O) 




10.15 < 10.35 < 10.46 


0.31 


0.0308 


0.32 


10.73 < 10.81 < 11.00 


0.27 


0.0450 


0.45 


9.27 < 9.34 < 9.51 


0.24 


0.0140 


0.04 




Geneva-Bessell(/ = 


0) 


10.25 < 10.27 < 10.34 


0.09 


0.0330 


0.32 


10.54 < 10.78 < 10.94 


0.40 


0.0400 


0.46 


9.18 < 9.39 < 9.51 


0.33 


0.0160 


0.04 


S3 


Padova-Bessell(/ = 


0) 


10.26 < 10.28 < 10.32 


0.08 


0.0330 


0.32 


10.50 < 10.74 < 10.90 


0.40 


0.0400 


0.45 


9.18 < 9.40 < 9.50 


0.32 


0.0160 


0.04 


-3 


SFR 




NGC2362 








A UI1 








V> Pot 

n Jr er 








Model (Vega cal) 




dm 


Adm 


Unc 


U(B-V) 


dm 


Adm 


Unc 


U(B-V) 


dm 


Adm 


Unc 


U(B-V) 




Geneva(O) 




10.48 < 10.58 < 10.69 


0.21 


0.0150 


0.12 


7.77 < 7.87 < 8.03 


0.26 


0.022 


0.12 


11.77 < 11.78 < 11.85 


0.08 


0.0250 


0.57 


-2 


V^cllc Vcl- -Dcooclll U J 




10.51 < 10.67 < 10.70 


0.19 


n ni 

U.UIOU 


n in 

U. -LU 


7.89 < 7.98 < 8.16 


0.27 




nil 

U. J. J. 


11.77 < 11.78 < 11.84 


0.07 


u.uzuu 


U. <J^± 




Padova-Bessell(O) 




10.47 < 10.66 < 10.77 


0.30 


0.0120 


0.10 


7.81 < 7.98 < 8.15 


0.34 


0.026 


0.11 


11.63 < 11.65 < 11.74 


0.11 


0.0310 


0.52 




Geneva-Bessell(/ = 


o) 


10.49 < 10.57 < 10.70 


0.21 


0.0120 


0.10 


7.86 < 7.96 < 8.13 


0.27 


0.024 


0.11 


11.76 < 11.78 < 11.80 


0.04 


0.0270 


0.54 


o 


Padova-Bessell(/ = 


o) 


10.47 < 10.55 < 10.72 


0.25 


0.0120 


0.10 


7.79 < 7.92 < 7.99 


0.20 


0.025 


0.11 


11.67 < 11.74 < 11.75 


0.08 


0.0300 


0.52 




SFR 






X Per 






NGC1960 






NGC2547 








Model (Vega cal) 




dm 


Adm 


Unc 


E (B-V) 


dm 


Adm 


Unc 


E (B-V) 


dm 


Adm 


Unc 


E (B-V) 




Geneva(O) 




11.82 < 11.82 < 11.85 


0.03 


0.010 


0.52 


10.21 < 10.22 < 10.30 


0.09 


0.0170 


0.22 


7.85 < 7.93 < 7.97 


0.12 


0.012 


0.053 




Geneva-Bessell(O) 




11.79 < 11.82 < 11.88 


0.09 


0.005 


0.50 


10.27 < 10.35 < 10.46 


0.19 


0.0164 


0.20 


7.98 < 8.05 < 8.09 


0.11 


0.018 


0.038 




Padova-Bessell(O) 




11.77 < 11.79 < 11.86 


0.09 


0.012 


0.50 


10.17 < 10.29 < 10.43 


0.26 


0.0164 


0.20 


7.97 < 8.04 < 8.06 


0.09 


0.018 


0.038 




Geneva-Bessell(/ = 


0) 


11.80 < 11.82 < 11.85 


0.05 


0.011 


0.50 


10.21 < 10.23 < 10.31 


0.10 


0.0164 


0.20 


7.92 < 8.03 < 8.07 


0.15 


0.020 


0.034 




Padova-Bessell(/ = 


0) 


11.79 < 11.80 < 11.82 


0.03 


0.010 


0.50 


10.20 < 10.31 < 10.40 


0.20 


0.0164 


0.20 


7.92 < 7.99 < 8.05 


0.13 


0.020 


0.034 





Table 6. The distance moduli and E(B — V) values , derived using U — B B — V fitting, for each SFR using the different models and calibrations to Vega colours of zero or Vega 
B — V = -0.002 and U — B = -0.004. Uncertainties required to achieve 0.40 < P r (r 2 ) < 0.60 are also shown. 



Revised Q-method 



SFR 




NGC6530 






NGC2244 








A Ori 






Model (Vega cal) 




dm 


Adm 


Unc 


E (B-V) 


dm 


Adm 


Unc 


E (B-V) 


dm 


Adm 


Unc 


E (B-V) 


Geneva-Bessell(O) 




10.49 < 10.50 < 10.60 


0.11 


0.011 


0.33 


10.66 < 10.77 < 10.81 


0.15 


0.01 


0.44 


7.99 < 8.01 < 8.12 


0.13 


0.005 


0.10 


Padova-Bessell(O) 




10.47 < 10.48 < 10.58 


0.11 


0.011 


0.35 


10.68 < 10.76 < 10.84 


0.16 


0.013 


0.45 


7.93 < 8.01 < 8.09 


0.16 


0.007 


0.10 


Geneva-Bessell(/ = 


0) 


10.36 < 10.41 < 10.49 


0.13 


0.011 


0.34 


10.67 < 10.76 < 10.89 


0.22 


0.014 


0.42 


7.95 < 7.97 < 8.11 


0.16 


0.007 


0.17 


Padova-Bessell(/ = 


0) 


10.48 < 10.52 < 10.57 


0.09 


0.009 


0.34 


10.70 < 10.81 < 10.92 


0.22 


0.013 


0.43 


7.94 < 7.97 < 8.08 


0.14 


0.007 


0.12 



Table 7. The distance moduli and mean E(B — V) values, after application of Q, for each SFR using the different models, and calibrations to Vega colours of zero or Vega B — V = —0.002 
and U — B = —0.004. Uncertainties required to achieve 0.40 < P r (r 2 ) < 0.60 are also shown. 
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Figure 8. h Per distance fit. 

These results show that any model dependency in the 
extinctions and distances derived is small. A final statistical 
justification can be found from the results in Section [6] Here 
the systematic uncertainties required to achieve a r 2 of ap- 
proximately one do not significantly favour any particular 
model. 



7 RESULTS 

To simplify our discussion of the implications of our derived 
distances, we have adopted the results from the Geneva- 
Bessell isochrones. We display the resulting fits as Figures 
EJ |3] and ISIUTI In Table [8] we provide a comparison of the 
adopted distances wi th 68% confidence i ntervals, with the 
distances assumed in Mavne etal . (2007). The best fitting 
E(B — V) is also provided. 



7.1 Notes on results 

We would have liked to include the sub-group CepOB3b 
in this paper, but after application of the Q-method the 
resulting colour range for the stars available for fitting was 
prohibitively l ow for distance fitting. Literature de rivations 
of extinction ([Garrison! Il970l : iBlaauw et al.l [l959) rely on 
intrinsic colours derived from other isochrones so cannot be 
used to fit with the Geneva-Bessell isochrones. 



8 IMPLICATIONS 

We have now derived a self-consistent set of distances (and 
extinctions) to, in general, a higher precision than that ex- 
isting in the literature, with statistically meaningful uncer- 
tainties for these distances. We now discuss some of the key 
implications of both the individual distances and of the en- 
tire dataset. 



8.1 Individual Distances 

Of the SFRs studied in this work distances derivations for 
eight are of particular note. Here we have converted the 
distance moduli to a distance to allow more obvious com- 
parisons. 
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Figure 9. h Per E(B-V) fit. Here a systematic shift in E(B-V) 
as a function of colour (U — B) is evident, as in Figure [3] This 
is probably due to a differe nce in the phot ometric systems of 
iBessell et all |l998) and Slesnick et all (j2002T ) (see Section [5^1). 
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Figure 10. NGC1960 distance fit. 
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Figure 11. NGC1960 E(B - V) fit. 
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SFR 


iviavne et ai. (zuin ) 
dm 


Range 


This work 

dm 


range 


A dm 


A range 


E(B - V) 


the ONC 


8.01 < 8.38 < 8.75 


0.76 


7.91 < 7.96 < 8.03 





12 


-0.42 


-0.64 


w 0.40( 1 ) 


NGC6530 


10.48( 4 ) 


» 0.40( 4 ) 


10.15 < 10.34 < 10.44( 2 ) 





29 


-0.14 


-0.11 


0.32 


NGC6530 


10.48( 4 ) 


w 0.40( 4 ) 


10.49 < 10.50 < 10.60( 5 ) 





11 


+0.02 


-0.29 


0.33( 5 ) 


NGC2244 


10.55 < 10.72 < 10.87( 3 ) 


0.33 


10.68 < 10.89 < 10.94( 2 ) 





26 


+0.17 


-0.07 


0.46 


NGC2244 


10.55 < 10.72 < 10.87( 3 ) 


0.33 


10.66 < 10.77 < 10.8l( 5 ) 





15 


+0.05 


-0.18 


0.44( 5 ) 


NGC2264 


9.6( 6 ) 


« 0.40( 4 ) 


9.26 < 9.37 < 9.52 





26 


-0.23 


-0.14 


0.04 


NGC2362 


10.84 < 10.87 < 10.90 


0.06 


10.51 < 10.67 < 10.70 





19 


-0.20 


+0.13 


0.10 


A Ori 


7.73 < 7.90 < 8.07 


0.34 


7.99 < 8.01 < 8.12( 5 ) 





13 


+0.11 


-0.21 


o.io( 5 ) 


A Ori 


7.73 < 7.90 < 8.07 


0.34 


7.89 < 7.98 < 8.16(2) 





27 


+0.08 


-0.07 


0.11 


a Ori 


7.41 < 7.80 < 8.19 


0.78 


7.84 < 7.94 < 8.10 





26 


+0.14 


-0.52 


0.06( ? ) 


X Per 


11.69 < 11.70 < 11.71 


0.20 


11.79 < 11.82 < 11.88 





09 


+0.12 


-0.11 


0.50 


h Per 


11.69 < 11.70 < 11.7l( 8 ) 


0.20 


11.77 < 11.78 < 11.84 





07 


+0.08 


-0.13 


0.54 


NGC1960 


10.40 < 10.60 < 10.80( 9 ) 


0.40 


10.27 < 10.35 < 10.46 





19 


-0.25 


-0.21 


0.20 


NGC2547 


7.92 < 8.18 < 8.47( 10 ) 


0.54 


7.98 < 8.05 < 8.09 





11 


-0.13 


-0.43 


0.038 



Table 8. Distance moduli and range assumed in lMavne et al.l (|2007l ) and those derived in this work. Notes are as follows. (1) Individual 
extinctions fro m T P ff, value given is approximat e mean. (2) Distance d erived using mean extinction derivation from r 2 fitting. (3) 
Not included in lMavne et al.l ([2007), distance from lHensberge et al.l ((2000). (4) Typical uncertainties assumed as none were provided in 
literature source. (5) Resulting distance after application of extinctions derived using Q-method (see Section |5?2T) . A mean extinction is 
quoted. (6) Typical uncertainties assumed as literature distance is the mean of many va lues. (7) E(B - V) fromBrown et al. | (ll994T). (8) 
Assumed to be at the same distance as x. P er in 

[Mavne et all (l2007f). (9) Not included i n lMavne et al.l (120071 ). distance from lPark fe Sund 
(2002). (10) Not included in lMavne et al.l ([20071 ). distance frorn lRobichon et al.l ([19991 ) . 
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Figure 16. NGC2547 distance fit. 




Figure 19. a Orionis distance fit. 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 0.3 



Figure 17. NGC2547 E(B - V) fit. 
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Figure 18. The ONC distance fit. Figure shown in intrinsic 
colour and extinction corrected magnitude (see Section [5j. 
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Figure 20. NGC6530 distance fit using Q-method extinctions 
(see Section |5.2|) . Figure shown in intrinsic colour and extinction 
corrected magnitude (see Section [5j. 
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Figure 21. A Ori distance fit using the Q-method extinctions 
(see Section |5.2|) . Figure shown in intrinsic colour and extinction 
corrected magnitude (see Section [5]). 
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Figure 22. NGC2244 distance fit using Q-method extinctions 
(see Section |5.2|) . Figure shown in intrinsic colour and extinction 
corrected magnitude (see Section [5}. 
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Figure 23. The ONC with a Geneva-Bess e ll 1 M yr MS isochrone 
and the Pre-MS isochrones of ISiess et all (l2000h for 1, 3 and 10 
Myrs. Stars appear to lie at the turn-on for an age of 10 Myr. 



8.1.1 The ONC 

We have increased the precision of the distance estimate 
for the ONC by a factor of 7 compared to that used in 
iMavne et all (|2007h . This new distance is also closer than 
the prev iously accepted res ult from the maser measure- 
ments of lGenzel etaP (|l98ll ). by 0.42 mag. Conversely this 
new distance, 391+g 2 pc agrees superbly with several re - 
cent derivations in the literature. Firstly, IJefTriesI (|2007ah 
finds a distance of 392 ± 32 pc from the rotational prop- 
erties of low-mass pre-MS stars (after removing accreting 
objects). Secondly, a parallactic distance of 389^21 P c from 
very long baseline arr ay observ ations has b e en fo und by 
ISandstrom et al.1 (|2QQ7h . Lastly, iKraus et all (|2007l ) find a 
distance of 434 ± 12 or 387 ±11 pc by modeling the orbit of 
the ^Ori C binary system. They adopt 434 pc as t he likely 
result a fter comparison to the distance obtained by I Jeffries! 
(2007a) for all objects including those showing evidence of 
accretion. Clearly, our result favours the solution yielding 
387 pc. A convenient round number which agrees with the 
majority of the recent derivations is 400 pc. 

This closer distance (391+g 2 pc compared to 480 =b 80 
pc) has important implications for the stellar population of 
the ONC. It means the pre-MS population lies 0.42 mags 
fainter in absolute magnitude in the CMD. This will force 
the isochronal age of stars older, but perhaps more impor- 
tantly increase their spread in isochronal age derived from a 
CMD (see the discussion in Pall a et al. 2005). This is due to 
the bunching of older isochrones towards the zero- age- main- 
sequence (ZAMS). In fact there is also evidence of a spread 
in the CMD from the MS members we have used to derive 
a distance. As can be seen in Figure 1231 MS stars exist at 
a position in the CMD suggesting isochronal ages of up to 
10 Myrs, which is at variance with the median pre-MS age 
of ~ 2 Myr (after allowing for the revised distance). This is 
discussed further in Section [8.2.21 

8.1.2 a Orionis 

We have improved t he precision of the distance estimate 
used for a Orion is in IMavne et al.l (|2007h by a factor of 2. 
ICaballerol (|2007l ) derives a distance based on TYCHO pho- 



tometry of 360lgQ pc. Although we also use the TYCHO 
photometry, w e use updated photometric conversions from 
iBesselll (|2000). We derive a distance of 389^4 

p c, a value 

more precise than, but in agreement with that of ICaballerol 
(|2007h . 



8.1.3 NGC2547 

For NGC2547 we have increas ed the pre cision in distance 
from that adopted in IMavne et all (|2007h by a factor of 5. 
In this work we derive a distance of 407^ 3 pc, which com- 
pares fav ourably with the HI PPARCO S result of 433^g 
c from Robichon et al. (1999). However, iNavlor &; Jeffries! 



20061 ) use pre-MS isochrones to obtain a distance of 361_ R 
c. T he age and distance derivation in INavlor fe J effries 
20061 ) is consistent with the Li depletion boundary, which 
is also based on pre-MS models. Therefore we conclude the 
difference in distance is attributable to a model-dependent 
difference between the MS and pre-MS models. Interestingly 
we als o obtain a d ifferent reddening, E(B — V) = 0.038 to 
th at oflClarial (1 19821) < E(B - V) = 0.06 ± 0.02 which is used 
in INavlor &; Jeffries! ([2006) . However, this does not have a 
significant impact on their distance, and so does not explain 
the discrepancy between the MS and pre-MS distances. 



8.1.4 NGC2244 

The distance to NGC2244 derived in this work is 1 425±?q 
pc, using individual extinctions from the revised Q-method 
(see Section I5.2|) . This compares well to th e literature result 
of 1 39 ±100 pc from eclipsing binaries of IHensberge et ahl 
(|2000l ). Previous MS is ochrone fitting s tudies placed this 
SFR at 1 667+ j 2 * pc (|Perez et all Il987l ). and 1660 pc 
(jPark fe Sung! i 2002. we use the same dat a as this study). 
Our result confirms the closer distance of IHensberge et ahl 
(2000) and is marginally consistent with the studies yielding 
a greater distance. The new distance will move the cluster 
pre-MS fainter or older by ~ 0.3 mags. 
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8.1.5 NGC2362 

Our derived distance for NGC2362 , 1 36llo7 p c is c loser 
and less prec i se th an that used in iMavne et al.l (|2QQ7h . In 
Mav ne et we adopte d a distance of 1493^9 pc, 

an incredibly precise value from lBalona &; Lanevl ((1996), de- 
rived using a form of MS fitting to narrow band photometry. 
For the sake of a consistent relative experiment we adopt our 
new distance for subsequent analysis. 



8.1.6 h and x Per 

The distance we have derived for \ P er °f 2 312+32 pc 
agrees well with the recent literature result of 2 344+53 pc of 
ISlesnick et al.l (|2002h , as stated in Section 14.11 They employ 
spectroscopy as well as the Q-method (with what appears 
to be an updated MS line, but canonical reddening vectors) 
combined with MS isochrone fitting to derive this distance. 
Interestingly they also derive a spectroscopic parallax dis- 
tance of 3 162 pc. This method involves dereddening stars 
with known spectral types onto an intrinsic MS; they con- 
clud e that recalibration o f these intrinsic colours is required. 
Also Sle snick et al. 1 (|20Q2h find h and x P er to be at approxi- 
mately the sa me distance, agreein g with the majority of the 
literature fe.g iKeller et al.ll200lf ). Our distance derivation 
for h and x P er are a l so consistent with both these clusters 
being at the same distance. 



8.2 Global issues 

In this section we discuss general implications related to the 
dataset as a whole. 



8.2.1 Metallicity 

Ther e is little work on the metallicity of SFRs. Ijames et al.1 
(2006) show that the metallicity is solar or very slightly sub- 
solar for the majority of stars in the Lupus, Chamaeleon 
and CrA SFRs. Conversely recent eclipsing binary re- 
sults suggest approximately ha lf-solar metallicity for h 
and y Per and Collinder 2 28 (jSouthworth et al.l l200l al: 
ISouthworth &; Clausen 2007), and solar metallicity for 
NGC6871 (|Southworth et al.ll2004bL Z = 0.02). If compo- 
sitions do indeed vary as is suggested from the above results 
then distances derived to these SFRs must be re-derived 
after a comprehensive composition survey. As discussed in 
Section 14.2.31 adopting a half-solar composition for x P er 
(for example) results in a fall of the derived distance mod- 
ulus by w 0.5 mags. A similar result would apply for h Per. 
This would have a severe impact on any age ladder, and thus 
on any conclusions about secular evolution, such as disc life- 
times. 



8.2.2 Age spreads and the R-C gap overlap 

As can be seen in Table [T]t he faint red edge of the MS defined 
using pre-MS isochrones and the edge of the apparent MS 
in the CMD do not agree for the younger SFRs. This means 
that there are MS stars which lie fainter and redder than 
the base of the MS predicted by theory. The models predict 
that, for a coeval population, the brightest stars still on 
the pre-MS (at the red edge of the R-C gap) are similar 



in magnitude to the faintest stars on the MS (which are at 
the blue edge of the R-C gap). This extension of the MS 
to magnitudes fainter than the head of the pre-MS we term 
the R-C gap overlap. The R-C gap overlap is most apparent 
in the younger SFRs. One of the best examples is shown in 
Figure 1231 the faintest stars on the MS in the ONC must 
be at least 10 Myrs old to have reached the MS, whilst the 
median age of the pre-MS is around 2 Myrs. An apparent 
age spread can also be seen for NGC2264 in Figure [25] for 
both the MS and pre-MS populations. 

This apparent extension (along a MS line) below the 
turn-on for almost all SFRs can be interpreted as an age 
spread within the SFR. This suggests that stars have evolved 
across the R-C gap before theory would predict for a coeval 
population and must therefore be older. This would explain 
why the R-C gap overlap is less apparent in the older SFRs. 
As the age of an isochrone increases they become fainter 
and move towards the ZAMS. Lower mass stars on the pre- 
MS contract more slowly as they age. Therefore, the same 
difference in age for an older population produces a smaller 
change in V than for a younger population. Thus the R-C 
gap overlap for older SFRs is harder to detect. 

Despite the R-C gap overlap being harder to detect in 
old SFRs, an isochronal age spread can still be seen in the 
HM stars. As shown for x P er m Figure [24] stars again lie 
below the turn-on and above the apparent turn-off for the 
best fitting or literature age of 13 Myrs. If one assumes the 
age is incorrect and increases it to match the turn-on, the 
turn-off will move fainter and exacerbate the problem for the 
brighter stars. Photometric variability and errors, and bina- 
rity have been shown not to completely account for these 
isochr onal age spreads for the pre-MS in iBurningham et al.l 
(|2005h and cannot account for the R-C gap overlap. 

Hypothesising such an age spread supports other results 
from isochrone modeling. It is well known that individual 
ages derived using pre-MS isochrone fitting also show an 
age s pread (s ee for example Palla et aT] [2005; Sles nick et al.l 
120041 ). Also, IJeffriesI (|2007bh finds a direct spread in the 
radii (and hence by implication age) of the PMS stars in 
the ONC, at a given effective temperature, a method free 
from isochrone theory. 

However it is dangerous to interpret the R-C gap over- 
lap, or the spreads in stella r radii and i sochro nal age, as real 
age spreads. As shown in iTout et"aL (1999) accretion can 
act to force a star bluer and temporally older. This shift in 
position within a CMD means an isochronal age spread de- 
rived from pre-MS fitting does not necessarily imply a real 
underlying age sp read of the same magnitude. Additionally 
ISiess et"aU |l999) find that the evolution of an accreting 
star is accelerated, with the star having a smaller radius 
and therefore lower luminosity than a non-accreting coeval 
counterpart. This result shows that age spreads derived from 
spreads in radii or from the R-C gap overlap again do not 
necessarily imply a real spread in age. 

Whether star formation is rapid (~ 1 Myr) or slow 

(~ 5—10 My rs) is c urrently an active debate (e.g 

jKrumholz &; Tanl 120071 ; iBallesteros-Paredes Sz Hartmannl 
[2007|). Not withstanding different star formation rates or 
several episodes of star formation, the age spread of the 
bulk of the population within an SFR is an approximate 
measure for the local star formation time. Apparent age 
spreads within a CMD are often used to support the model 
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of slow star f ormation (Pal la et al.ll2005l : iBurningham et all 
2005; Jeffries 2007b|). In the rapid star formation model these 
spreads of apparent age are dominated by accretion effects, 
and the residual real age spread is small. If accretion does 
indeed act to scatter a star within the CMD and even ar- 
tificially accelerate its evolution or contraction, isochronal 
ages should not be used to represent real age spreads. More- 
over, isochronal ages (based on turn-ons, turn-offs or pre-MS 
fitting) for individual stars without a known accretion his- 
tory do not represent the true age of the star and therefore 
should not be used to support evolutionary theories. Indeed, 
it is even hard to argue that a median or mean age for a 
given SFR derived from isochrone fit ting has any real m ean- 
ing. Perhaps, foll owing the results in lSiess et al.l ((1999) and 
iTout et all (jl999h , if accretion causes a decrease in the star's 
radius, therefore increasing its isochronal age, the most ac- 
curate representative age for a given cluster is that of the 
youngest stars having the lowest accretion histories. How- 
ever, this would mean a dramatic change in the ages for 
most SFRs, for example the yo ungest stars in the ONC are 
w 0.1 Myrs (|Hillenbrandlll99/t ) and many older SFRs still 
contain active star formation and embedded objects. 

A useful indicator to quantify these perceived spreads 
may be the R-C gap overlap, where the minimum mass ob- 
ject to have developed a radiative core and to have joined 
the MS can be compared to the maximum mass stars still ex- 
isting on the convective pre-MS. This overlap then provides 
a precise diagnostic for the apparent age spreads. Whether 
these spreads show a real age spread or are indicative of the 
range of accretion histories present depends on the model 
adopted. 

Given the problems with SFR ages the best approach 
currently is to compare observations of two different SFRs. 
Either deriving an age order by assuming a similar range of 
accretion rates within each SFR, or by using the R-C gap 
overlap to derive approximate differences in the range of ac- 
cretion rates. An example of this can be seen in Figure 1251 
showing the absolute magnitude and intrinsic colour for the 
stars of the ONC and NGC2264 (taken from iMavne et all 
2007, and sources referenced therein) . The locus of the pre- 
MS in NGC2264 clearly lies slightly below that of the ONC, 
but the MS section is strikingly similar. In addition Figure 
[25l shows the MS for the SFRs extending below the pre- 
dicted turn-on. Moreover the brightest MS stars in both 
populations are at similar magnitudes i.e. the turn-off is in a 
similar position with stars lying on the apparent MS above 
the turn-off (as seen in Figure [2^]) . As the apparent MS and 
pre-MS in both clusters appear to extend to similar points, 
i.e. the R-C gap overlap is similar, suggesting a similar range 
of accretion histories or ages for each cluster. So, in conclu- 
sion, the MS and pre-MS sections imply a large isochronal 
age spread as seen in Figure [23] and deriving ages and age 
spreads from these sections of the sequence would lead to 
a similar result for the ONC and NGC2264. However, com- 
paring the sequences as a whole show that NGC2264 is more 
evolved than the ONC. 

8.2.3 Secular evolution and ages 

We now reconstruct the age ladder of IMavne et all (|2007h 
using the new distances and extinctions from Table [8] We 
have however changed the process slightly. Given that the 
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Figure 24. x P er with a Geneva-Bessell 13 Myr MS isochr one 
(dashed line) and the pre-MS isochrones of lSiess et al.l (l2000h for 
13 and 23 Myrs (bold lines). Stars appear to lie below the turn-on 
and above the turn-off for an age of 13 Myr . 



NGC2264 and the ONC 
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Figure 25. Absolute magnitude as a function of in trinsic colour 
for s tars in the ONC and NGC2264 (taken from IMavne et all 
2007, and sources referenced therein) adjusted using the extinc- 
tions and distances from this work. The blue dots are the ONC 
and the red dots are NGC2264. The locus of points for NGC2264 
can be seen to lie below the ONC, however the MS section is very 
similar in extent for each SFR, i.e . the turn-o n appears coinci- 
dent. The pre-MS isochrones from Siess et al. (2000) are shown 
for 1 and 3 Myrs. Most importantly in the MS section stars from 
both populations extend below the turn-ons whilst for both SFRs 
the peak of the MS (brightest stars) is similar i.e. the turn-offs do 
not appear significantly different, although this is based on only 
a few stars. 

new distances are generally more precise, we assume that 
the MS for each SFR will be approximately coincident in the 
CMD. Therefore, we plot the individual photometric points 
for this section of each sequence. Then for stars redward of 
the R-C gap we fit a spline through the median points. We 
retain the ZAMS subtracted space as a presentation tool 
(where the colour at each magnitude h as the correspondin g 
ZAMS colour subtracted from it as in IMavne et al.ll2007h . 
We are unable to inc lude NGC2244 as t he pre-MS of this 
SFR is not studied in Mavn e et al . (2007). We do however, 
include in this section four SFRs from Mavn e et al . (2007) 
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Figure 26. The MS stars (dots) and pre-MS empirical isochrones 
(continuous lines) of the ONC (red), h and x P er (black) 
and NGC2362 (green), with the pre-MS empirical isochrone of 
NGC2264 (blue) as a dashed line. The left panel is a CMD of 
absolute magnitude and intrinsic colour with the right panel is in 
ZAMS subtracted space. 
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Figure 27. The MS stars (dots) and pre-MS empirical isochrones 
(continuous lines) of NGC2264 (blue) and h and % Per (black), 
with the pre-MS empirical isochrones of A Ori (red) and a Ori 
(green) as dashed lines. The right panel shows the same in ZAMS 
subtracted space. The left panel is a CMD of absolute magnitude 
and intrinsic colour with the right panel is in ZAMS subtracted 
space. 



not studied in this work, where we adopt literature distances 
and extinctions. 

As the distance moduli for some SF Rs have changed 
from those adopted in Mavn e et al . (2007) the nominal ages 
may also have changed. The distance moduli which have 
changed by more than 0.2 mag which were discussed in 
iMavne et all (|200/t ) are the ONC, NGC2264, NGC2362, A 
Ori and a Ori. In addition the relative distances and extinc- 
tions for h and x P er have changed sli ghtly, therefore we 
have refitted these clusters. Following IMavne et al.l (|2007h 
we have adjusted h Per to the distance and extinction of x 
Per and combined the sequences prior to fitting them. We 
have included the entire sequence in a fit to h and x P er 
providing an e mpirical ZAMS blue -ward of the R-C gap. 

Following IMavne et al.l (|2007h we present fiducial em- 
pirical isochrones bounding the target empirical isochrone 
in a CMD of absolute magnitude and intrinsic colour and 
in ZAMS subtracted space. Figure [26] shows the pre-MS of 
NGC2264, after application of the distances in this work, 
lying only slightly below the ONC, the pre-MS of NGC2362 
is shown to lie below these sequences. A combined h and 
X Per empirical isochrone is also shown as a lower fiducial 
in Figures [261 and [271 with Figure [271 using NGC2264 as an 
upper fiducial. Figure [271 also shows the empirical isochrones 
for the pre-MS of both A Ori and a Ori, the former lying 
above NGC2264 and the latter below. From these two fig- 
ures we can create an age ladder (youngest to oldest) and 
assign nominal ages: the ONC (2 Myrs), A Ori , NGC2264 
and a Orionis (3 Myrs), NGC2362 (4-5 Myrs) and finally h 
and x P er (13 Myrs). 

We have repeated this method to obtain positions in 
an age ladder for each SFR for which we have enough 
data to create an empirical isochrone. This has resulted 
in the creat i on of a new group in addition to those of 
IMavne et al . (2007). This group, at a nominal age of 2 Myrs, 
contains the ONC and NGC6530, as it lies older than IC5146 
but marginally younger than NGC22 64. The result i ng ag e 
groups from this work and those from IMavne et al.l (|2007h , 



alongside literature estimates of the fraction of stars exhibit- 
ing infrared (IR) excess (i.e. candidates for an associated 
disc) , are shown in Table [9] It is important to note that we 
have only updated distances for some of the SFRs. 

The infrared excess fractions from Table [9] can be used 
to infer the presence of a disc, then by comparing the disc 
fraction across a range of SFRs one can ex amine the evolu- 
tion of these discs, as in Haisc h~et al.l (200l|). Figurel28lshows 
the logarithmic nominal ages for the SFRs from Table [9] and 
the inferred fraction of stars with discs. However the disc 
fractions inferred from the data in Table [9] come from dif- 
ferent mass ranges, dependent on the apparent magnitude 
range (and therefore distance) of the target SFR. Addition- 
al ly, the excess criteria us ed are different, with studies such 
as lHernandez et al.l (|2007T ) using the Spitzer Space Telescope 
IRAC and M I PS ca mera channels, whereas disc fractions in 
Haisc h et al . (2001) were calculated using JHKL excesses. 
Therefore, given the heterogeneous nature of these data one 
cannot draw strong conclusions from Figure [28] regarding 
the disc fraction as a function of age. 

Despite its limitations, Figure [28] is not consistent with 
a uniform decay, revealing further possib le evidence of envi- 
ronmental effects as suggested in lMavne et al.l (|2007h . As an 
example we examine the inferred disc fractions of three SFRs 
in the same age group (nominal age of 3 Myrs); NGC2264, 
A Orionis and a Orionis, with distances of dm = 9.37, 8.01 
and 7.94 respectively. The disc fra ctions adopted are (in 
the same order) 52 ± 10% (JHKL, lHaisch et all l200ll . for 
masses greater than 0.85Mp), w 25% (from IRAC data, 
iBarrado v Navascues et all 120071 . all discs in spectral range 
of M0-M6.5 or approximate mass range of 0.1 — 0.8Mq 
using pre-MS isochrones), a nd 31.1 =b 3.8% (from IRAC 
data, Hernand eTet all [2007, all stars in the approximate 
mass range 0.1 — I.OMq). In the case of these three SFRs 
NGC2264 has an inconsistent disc fraction, it is much higher 
than that of the other two SFRs. These disc fractions are 
taken from differing mass ranges and the SFRs are at differ- 
ent distances which could lead to sensitivity problems in the 
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SFR 


This work 


Nominal a 
Mavne 


ge 

et al. (2007) 


Fraction of stars with IR excess 


IC5146 


- 




1 


- 


NGC6530 


2 




1 


44% (7) 


the ONC 


2 




1 


80 ± 5%W 


A Ori 


3 




3 


« 25% « 14% « S1%( 6 ) 


CepOB3b 






3 




NGC2264 


3 




3 


52 ± 10%^) 


a Ori 


3 




4-5 


31.1 ± 3.8%, 26.6 ± 2.8% and 33.9 ± 3.1%( 2 ) 


NGC2362 


4-5 




3 


12±4%( 1 ),7±2%( 3 ) 


IC348 






4-5 


65dz8%( 1 ) 


NGC7160 






10 


« 20% W 


h and x P er 


13 




13 


2 - 3%( 8 ) 


NGC1960 


20 




20 


3 ± 3%^) 


NGC2547 


40 




38 


« 7%( 5 ) 



Table 9. The new nominal ages of SFRs from this work and those fro m I Mavne et al.l (2007) are shown in addition to the fraction 
of stars with IR excesses. NG C2244 from this w ork was not studied in Mavn e et al.l (|2007h so is omitted. IC5146, CepOB3b, IC348 
and NGC7 160 are studied in | Mavne et al.l (Eo07f) and includ ed here by adopting literature distances and extinctions. The notes are as 
follows. (1) lHaisch et al ] (l2Q0lh . (2) iHernandez et al.l (|2007l ). First value for TTS stars (approximate mass range 1 — O.IM^). Se cond 
and third values for entire sam ple, first for t hick discs and second for thick and evolved discs. (3) iDahm fc Hillenbrand! (|2007h . (4) 
ISicilia-Aguilar et al.l (I2005I V (5) lYoung et al.l (I2004T ). (6) iBarrado v Navascues et al.l (I2007I V Firstl value derived from the IRAC CCD. 
Second and third from spectral-energy di stribution (SEP) fitting, thi ck discs then thin and thick discs combined. All values in the 
approximate mass range 0.1 — I.OMq. (7) IPrisinzano et al.l (|2007l ). (8) ICurrie et al.l (|2007l ). for h and x P er stars fainter than J=13.5, 
stars brighter than this have 0% disc fraction. It is important to note that each of these studies is over a different mass range (due to 
different distances and apparent magnitude ranges) and excess candidates were selected in a heterogeneous fashion. 
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Figure 28. Figure showing the log(nominal age) of each SFR 
with the IR excess from various literature sources. The inset panel 
shows the region 0-5 Myrs enlarged. 



L band as suggested for NGC2362 bv lLvo et all f|2QQ3h . How- 
ever, in this case the further distance to NGC2264 would re- 
sult in fewer L band detections and the lower mass limit be- 
ing higher is also likely to decrease the detected disc fraction. 
Therefore, it is likely that a consistent experiment would in- 
crease the discrepancy between NGC2264 and the two SFRs 
with lower disc fractions. However, even ignoring the partic- 
ular case at a nominal age of 3 Myrs it is clear that the these 
data do not necessarily imply a smooth decline in disc frac- 
tion with age, suggesting other, presumably environmental 
factors may affect disc lifetimes. 



(i) We have derived a self-consistent set of distances of 
generally higher precision than previously available for a set 
of SFRs in the age range 1-40 Myrs. We have also derived 
distances using several other models and calibrations (see 
Section [6]). 

(ii) In addition to these new distances and reddenings 
(or extinctions) w e have reconstructed the age ladder of 
iMavne et ail (2007) assigning new nominal ages, as shown in 
Table [9] To enable the reader to add other SFRs to this lad- 
der the pre- MS splines are freely available from the cluster 
collaboration home pag^E an d the CDS archive. 

(iii) We have shown that metallicity information is now 
vital for accurate relative distances to SFRs. This is espe- 
cially true if one is attempting to characterise evolutionary 
indicators such as disc fractions as a function of age, or try- 
ing to uncover environmental effects (such as the effect of 
ionising winds from massive stars on planet formation from 
discs). 

(iv) We have discussed that the overlapping region of the 
R-C gap could be used to derive spreads in isochronal age in 
a CMD, i.e. the apparent age difference between the maxi- 
mum mass star still on the convective pre-MS and the mini- 
mum mass star which has reached the MS. If star formation 
is slow and isochronal ages of individual stars are reliable 
this would provide a direct measurement of the age spreads 
present in SFRs. If star formation is rapid the R-C gap over- 
lap region reveals the underlying spread in accretion histo- 
ries within an SFR. This is important as for rapid star for- 
mation, if an accretion history is unknown isochronal ages 
derived from a position in a CMD do not represent the true 
age of a star. Indeed it is therefore likely that if a rapid star 



http://www.astro.ex.ac.uk/people/timn/Catalogues/description.html 
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formation model is accurate median or mean ages drawn 
from a population are also invalid. A more useful approach 
may be to compare SFRs using age ladder arguments or 
perhaps to use the age of the youngest stars which have the 
lowest accretion history. 

(v) We have shown further evidence for non-uniform de- 
cay of discs in SFRs, although new comparisons must be 
made using consistent disc fraction indicators and mass 
ranges. 
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